o 
o 

(N 



Oh 
O 



The Dynamical Mixing of Light and 
Pseudoscalar Fields 



Sudeep Das a , Pankaj Jain 6 , John P. Ralston c and Rajib Saha 5 
"Department of Astrophysical Sciences 
Princeton University 
New Jersey 08544, USA 
^Physics Department, IIT, Kanpur - 208016, India 
\£> ; c Department of Physics & Astronomy 

^ ! University of Kansas 

Lawrence, KS-66045, USA 

O ; February 2, 2008 



Abstract 



We solve the general problem of mixing of electromagnetic and 
scalar or pseudoscalar fields coupled by axion-type interactions Ci nt = 
^ . g^cj) e^apF^ F a ^ . The problem depends on several dimensionful scales, 

including the magnitude and direction of background magnetic field, 
the pseudoscalar mass, plasma frequency, propagation frequency, wave 
number, and finally the pseudoscalar coupling. We apply the results 
to the first consistent calculations of the mixing of light propagating 
in a background magnetic field of varying direction, which shows a 
great variety of fascinating resonant and polarization effects. 

For about 20 years the mixing of light and pseudoscalar fields in 
propagation has been studied with fascination P-JZj. The subject 
generated renewed attention in the context of cosmological observ- 
ables that can probe exceedingly small couplings [HI 03 E3 IH] ■ One- 
recent approach proposes that the dimming of supernova light might 
be explained by transition of light into unobserved pseudoscalar, or 
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"axion," modes although this effect might be limited by obser- 
vations of radio galaxies |13| . It has also been pointed out that pseu- 
doscalar field can generate magnetic fields due to their coupling with 
photons |14| . Polarization observables are even more sensitive than 
intensity: for coupling constants many orders of magnitude too small 
to cause dimming, the cumulative evolution of phase shifts can gener- 
ate phenomena clearly violating the Maxwell equations in plasmas |15j. 
Several laboratory experiments have also sought the spontaneous res- 
onant conversion of dark matter axions to photons, and explored the 
possibilities of conversion in lab-made magnetic fields. 

There is a well-established theoretical technology of mixing light 
with a background magnetic field transverse to propagation. Yet de- 
spite long study, we know of no complete solution to the mixing prob- 
lem depending on every possible variable. And there is no wonder, 
as there are many dimensionful scales, including the magnitude and 
direction of background magnetic field, the pseudoscalar mass, plasma 
frequency, propagation frequency, wave number, and finally the pseu- 
doscalar coupling. By approaching the problem with new methods 
here, we will be able to survey various limits used in the literature 
and also present a convincing resolution of the dynamics in a slowly 
varying background field of arbitrary direction. 

The basic Lagrangian assumes a pseudoscalar 1 field 4> coupled to 
the electromagnetic field strength F^ u by the action 



We include a coupling to a current for completeness. For the pur- 
poses of linear propagation the potential V(cj)) can be ignored as a 
small perturbation, and the metric g replaced by a given background 
form. Certain non-local plasma effects, described by the plasma fre- 
quency, Faraday rotation, etc., may also need to be incorporated. 

By translational symmetry, certain eigenmodes will evolve like e tkiZ 
in propagation over a distance z, where hi are wave numbers to be 
determined. This is simple and obvious. Yet one might claim the 
opposite that k should be fixed, while frequency uj remains to be 
determined, as so common in quantum mechanics and neutrino os- 
cillations. Indeed some literature solves for uj eigenvalues without 

We may let cf> also be a scalar field given parity violation. 
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discussion. Physics is local, and by Huygen's principle, i.e. the use of 
causal Green functions, a source with known time dependence er lu)t 
must over distance develop its own wave numbers fcj to propagate 
on-shell. There are no boundary conditions of fixed k from which to 
calculate frequency, so it is not negotiable that fcj are to be solved, 
just as in careful work on neutrino oscillations ^3 El- We also give 
extra attention to maintaining gauge invariance, which we have not 
seen before. The physics turns out to be surprisingly intricate. 

We apply the revised propagation equations to the interesting 
problem of light traveling in a background magnetic field of varying 
direction. For the parameter values of axion masses, magnetic fields 
and couplings commonly assumed, the magnitude of new changes is 
often non-negligible. This fact aside, the results themselves are fas- 
cinating, and full of remarkable complexity and structure, somewhat 
like a generalized version of the resonant propagation of neutrinos. We 
think this is very interesting: The possible existence of axions can be 
probed in polarization observables for parameters ranges far smaller 
than will cause a dimming of light by direct conversion. Although 
axion-related dimming is given some credence it is usually assumed 
there are no exotic polarization effects to be observed. We find that 
the absence of exotic polarization effects would be able to rule out the 
light-dimming hypothesis. Confrontation with data on polarization, 
of course, needs a detailed study of many potential backgrounds to 
any signal, and would go beyond the scope of this paper. Our main 
task is simply to get the propagation equations resolved once and for 
all. 



1 Gauge Invariant Methods 
1.1 Equations for E and 4> 

To eliminate difficulties of gauge invariance we first obtain the non- 
covariant form of the Maxwell equations with no approximations JH] : 

V-E = 9< F<l>-{0+d)+p\ (3) 
VxB-^ = gjExV^-(B + B)^\+j; (5) 
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V •(*? + £?) = 0. 



(6) 



Here + B { = \e ijk Fi k and E { = F 0i are the usual magnetic and 
electric fields. Here B and B represent the magnetic field due to the 
background and due to the electromagnetic wave respectively. 

In anticipation we note that the revised "Gauss's Law" Eq. |3] cou- 
ples the longitudinal electric field to V(/>. This creates a qualitative 
change compared to light in free space, where the longitudinal mode 
does not normally propagate. If (ft propagates we now have a propa- 
gating longitudinal light field. If there is a plasma, then the ordinary 
Gauss's Law becomes V • eE = pf ree , where e is the dielectric constant 
(or "permitivity ") . Since e = e{to) is not local in the time domain, we 
will incorporate it below in the Fourier-transformed equations. 

The pseudoscalar field's equation of motion is 

^ - V 2 + m%4> = - 94> E ■ (B + B) (7) 

Gauge invariance is explicit, and one can check current conservation 
directly, 

Assume B solves the zeroeth order Maxwell equations with no (ft 
background. The linearized equations for E/c « B, B « B are 

V-E = g^VcP-B + p; (8) 
dB 

VxE + — = 0; (9) 

- dE -?d(ft , . 

V-B = 0. (11) 

Proceed to get a wave equation for E by taking the curl of Fara- 
day's Law, 

VxVxB = -V 2 E + VV • E = — ^-V x B, 

ot 

and substituting Eqs. Replacing B ~ B gives 

" V 2 ^ + ^ = ~ ^V(V0 • B). (12) 
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In this equation the longitudinal part of E mixes with V</>. Take the 
transverse(sub-T) and longitudinal parts (sub-L) of the electric wave 
equation, for wave number k, with El = k ■ E: 

(k 2 -u; 2 )Et = -g^&rfa (13) 
{k 2 -u J 2 )E L = g <t> {k 2 -u J 2 )BLct>. (14) 

There clearly exists no gauge in which the longitudinal electric field 
decouples from the problem. If we limit the study to ■ B = 0, 
then Gauss' Law makes E transverse. Everything in the literature is 
perfectly consistent. 



1.2 Equations for D and 

Another method is needed when k ■ B ^ 0. Many linearized electro- 
magnetic theories can be encompassed by the equations: 

0; (15) 

0; (16) 

0; (17) 
0. (18) 

The purpose of the "archaic" representation via D is to have a field 
which is perfectly transverse. With D the transverse wave operator is 
greatly simplified: 

V x (V x D) -» -V 2 D. 

This effectively reduces the freedoms of the propagating gauge fields 
from 3 to 2: one would have to use 4-state mixing of 3 E components 
and one <j) if this were not arranged. 

Can we make D and H serve in Eqs. 151- 1111 and also include plasma 
effects? We find Eqs. 1151181 consistent with the definitions: 

D = eE-g^B; (19) 
H = B . (20) 

The asymmetry here comes from having a magnetic background. In 
our work we will assume the contribution to e due to the plasma 



V • D 

^ ^ dB 



V x H 



dD 
~ ~dt 
V -B 
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frequency u p , via 





1.3 Decoupling 

From Faraday's Law and the D equation we have 

1 V x (3 + g^Bcp) 

1 - (uj^/uj 2 ) 

VxVx(fl + ^)=[l--^)^. (21) 

Together with the <f) propagation from Eq. the equations have 
been simplified as much as generally possible: the coupled system 
of (j>, D x , D y , D z have one locally decoupled mode, no longitudinal 
mode, and are equivalent to two coupled pde's with no approximations 
other than linearization. 

We now drop terms of order X7B/B as negligible compared to other 
length scales, including the splitting of modes, setting up the usual 
adiabatic limit. We seek local plane wave solutions with V — ► ik. The 
component of D perpendicular to B decouples: 

{k 2 + u: 2 p -u 2 )D x B = 0. (22) 

The other transverse projection of the D wave equation becomes 

( k 2 + u/jj - to 2 )D ■ B T + k 2 g<t> BT^ = 0. (23) 

Notice that in using D the equation of motion, involving the curl 
of B, is not used: in fact it is satisfied as an identity. Conversely, 
when Faraday's Law is substituted into the E wave equation, then 
Faraday's Law is satisfied as an identity, and the equation of motion 
is solved (Eq. El)- By subtracting Eq. [221 from the (in principle) 
independent wave Eq. El for E at the compatible point, we obtain a 
nice consistency check. 

We turn to the coupled system: 

(k 2 +uj 2 -u 2 )D- B T + k 2 g^Brcl) = 0; (24) 
9 " BT D.B T +(k 2 + ml-u 2 + ^ #f, \ = 0. (25) 
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The system can be solved directly for the dispersion relation k 2 = 
k 2 {uj) by setting to zero the determinant of the corresponding matrix 
M, defined by 



However the eigenvalues k 2 needed are not on the diagonal. Moreover 
M is not symmetric, and non-symmetric matrices have eigenvectors 
which are not orthogonal. 

Much the same occurs in optics |19j . where the corresponding equa- 
tions for propagation with a tensor dielectric constant are: 



multiplication on the left by 6 (k) yields a symmetric eigenvalue equa- 
tion: 



The propagation eigenstates are obtained from the 2x2 matrix 5 T (k)e~ 1 5 T (k) 
in the sector transverse to k. This is considerably more subtle than 
(say) diagonalizing first, and simply taking a transverse part. 

This indicates that further transformations are needed for a useful 
solution. 




(k 2 S T (k)e- % -u 2 )D = 0; 
One seldom finds 5 T (k)e~ l to be symmetric. Yet since 



(26) 



(k)D = D, 



(k 2 6 T (k)e- 1 5 T (k)-u J 2 )D = 0. 



(27) 



1.4 Orthogonal Modes 

First, D x B decouples from (f> and propagates like ordinary light (in- 



cluding plasma frequency) with wave number ko = \ uo 2 — uj 2 




We made the rest of the transformation by inspection. Define 



D = 



D-B + gB T (p 



(28) 
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Now the propagation matrix is symmetric and eigenvalue k 2 lies on 
the diagonal: 



k 2 + ujp - u) 2 g<t>BTU \ ( D 



g4,B T u k 2 + m 2 ,- to 2 



with 



q 2 B 2 

™>l = ™l + i 1 2 1 2\ ■ ( 30 ) 

As a consequence propagation generates unitary rotations of (D, </>). 
Go to a new basis 

i i „ /cos(9 — sin0\ , . 

7? >= O \ip A >, O = [ . n a ( 31 ) 

V sin a cos a J 

The mixing angle diagonalizing propagation is 

tan20 = ^f^. (32) 



(29) 



The dispersion relations are 



where 



k\ = ^-IfmJ + ^-iVn*; (33) 

fcl = W 2_ 1 ( ^2 +(J 2 )+ 1^ (34) 



n* = 44^ W a + (mJ- W J) 2 . (35) 



By inspection of these results, the eigenvalues and mixing are just 
the same as solving the Bl = limit and making the replacement 
m 2 ^m 2 = m 2 +g 2 B 2 L /[l - (u 2 /u 2 )}. 

1.4.1 Plane Wave Simplification 

There are circumstances where neglecting X7B/B may be not possible. 
Then Eq. and Eq. [7| cannot be simplified further. However if 
the propagation can be reduced to plane wave modes with constant 
parameters, there is a simple way to understand the modes. 
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First solve the longitudinal mode using Gauss' Law: 

k-eE = g^U) ■ 
eE L = kg^U-B. (36) 

Here k = k/k is a non-local operator. Insert the solution where it 
appears in the propagation of 0, Eqn. 

(Lu 2 + k 2 + m 2 )0 = -g^E-B, 

-> -g,E T B T - 9 ^ (37) 

Observe that the effects on equations for <j) are the same as replacing 
m \ ~ * = m \ 9^l/ € - Meanwhile the transverse projection of 
the electric equation, Eq. 1121 also involves only Et and (p. Since 
this subsystem has decoupled, they must have modes which are linear 
combinations of (j) and Et- finally we recover the transformation to 
reveal that D = Ex in this limit. 



1.4.2 Non-Perturbative Effects, and A New Resonance 

The axion mass must be very small, and it often appears safe to take 
— > 0, uj p — > smoothly. This is the assumption of several phe- 
nomenological applications. Let us revisit that question in light of the 
mixing solution Eq. 1321 which reduces to 



- Bx 1 ^ 

tan 29 = ^, = tu p = 0. (38) 

9<f>Bc 

This is a very interesting result: the mixing is inversely proportional 
to the coupling constant, a typical non-perturbative effect. 

Can we take Eq. ESI seriously? In current cosmological propagation 
applications, the order of scales is lo 2 > m| > g"lB 2 , and Eq. 1351 does 
not apply. However in lab experiments we can control the vacuum to 
make lo p ~ to a very high degree. 

Restoring the full dependence on all variables, tan(29) of Eq. 1321 is 
compared to the traditional formula valid for Be = in Fig. ^ There 
is a type of resonance at lo = lo p , with width 

A ^ 

AlO ~ LO r , 
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Given current values of g^, the resonance may be too narrow to ob- 
serve. 

The resonant effect is very interesting conceptually. Qualitatively, 
at resonance it appears that the longitudinal mode of a plasma oscil- 
lation becomes very strongly mixed with the pseudoscalar field, de- 
pending on the difference of masses. As we mentioned earlier, (f> mixes 
with El: indeed due to the constraint of Gauss's Law, it is the same 
dynamical phenomenon as the longitudinal field. Let us estimate some 
magnitudes: when fully mixed, (f> ~ D, or 




As mentioned earlier El ~ 9 ' p - c . Together the relations predict 

Et eujp uj^ — uj 2 

Thus there is always a frequency for which we may observe the for- 
merly non-interacting pseudoscalar electromagnetically, and as a form 
of longitudinally polarized light : The El being observable and af- 
fecting instruments just as much as a longitudinal field in a plasma 
oscillation. Given sufficiently fine measurements the "invisible axion" 
could in principle be "visible." 

We hope to explore more deeply the potential laboratory repercus- 
sions of these phenomenon in another paper. Given that most current 
interest centers on cosmological propagation, we turn to studying the 
effects of a varying B field in the next Section. 

2 Three Mode Mixing: Varying B 

We next consider the adiabatic propagation of light through a back- 
ground magnetic field which varies slowly in direction. This problem 
has not been solved before. The results are far from trivial, and give 
substance to many cosmological applications assuming some "fluctu- 
ating" magnetic fields with typical coherence lengths. As we will show, 
the variety of physical phenomena one can observe is very great. In 
some limits, writing a transition probability and taking a statistical 
average may suffice, but in other limits the polarization effects are 
quite spectacular. The dynamical possibilities for the mixing of light 
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Figure 1: Behavior of tan(28) in the vicinity of to ~ uj p for m 2 ~ cu 2 . The 
dashed line shows the calculation for = 0. Parameter values have been 
rescaled to make the resonance visible: it would be exceedingly narrow given 
current beliefs for g^. 



actually exceed those for neutrino-mass mixing, which has been stud- 
ied for nearly 50 years and still appear inexhaustible. Indeed we were 
engaged in the general propagation problem by this very attractive 
method to probe the existence of light pseudoscalars in cosmology 

The physically observable density matrix p is given by 



< E\\Et > < E\\E* ± > 

< E±En > < E ± E* ± > 



(39) 



where < > denotes the statistical averages occurring in propagation 2 . 

Orient the z-axis along the direction of the wave. Let angle £ 
measure the direction of the background field relative to the x-axis: 

B T = B cos £(z)t + Bsin£(z)j (40) 

We fix the magnitude of the background magnetic field to identify 
effects arising due to varying magnetic field direction. A changing 



2 We decline to develop a 3 x 3 density matrix including the longitudinal mode, as 
unlikely to be observed in these circumstances 
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background magnetic field magnitude is easily included in the formal- 
ism. For the same reason we ignore the variation in plasma density 
along the path. The effects of varying plasma density for fixed field 
direction has been studied in detail elsewhere |15j . 
The wave equation can now be written as 

/ A x \ I uj 2 -gBucos^ 

(u; 2 +d 2 ) I A y - u? p -gBu; sin f 

\ <p J \—gBu> cos £ — gBujsmt; m\ 

(41) 

where A = E/u. We dropped g 2 B\ terms as negligible for inter- 
galactic propagation with typical parameters. With a slowly varying 
background and working in the adiabatic limit, we define transformed 
fields A' x , Ay and <p' such that 



(42) 




A x \ 




( cos (3 


— sin (3 




H 


sin (3 


cos (3 






V o 







The wave equation reduces to 
iu; 2 + d': 








^ 









uj 2 


-gBuj 


(?)- a 


(43) 





-gBu 









Here j3 = £ — n/2. The equation reduces to the case of two component 
mixing which can be solved along the lines discussed in Ref. |15j . 
Once we have obtained all the correlators between A' x and A' y , we can 
express the required correlators as 

< A* x {z)A x {z) > = sin 2 £(,z) < A'*{z)A' x {z) > +cos 2 £(z) < A'*{z)A' y {z) > 

+ cos£(.z)sin£(z) (< A'*(z)A' y (z) > + < A'*(z)A' x (z) >) 

< A* y (z)A y (z) > = cos 2 £(z) < A'*(z)A f x (z) > + sm 2 £(z) < A'*(z)A' y (z) > 

- cos£(z)sin£(z) (<A':(z)A'y(z)> + <A' y *(z)A' x (z)>) 

< A* x (z)A y (z) > = - cos t(z) sin az){<A' x *(z)A x (z)>-< A'; (z)A y (z)>) 

+ sm 2 az)<A' x *(z)A' y (z)> 

- cos 2 i{z)< A 1 * (z)A' x {z)> . (44) 

The correlators appearing on the right hand side of these equations 
can be calculated by using the results in Ref. [T5j . 



12 




0.1 1 10 100 1000 0.1 1 10 100 1000 




0.1 1 10 100 1000 



Figure 2: Normalized Stokes parameters (a) Q/I, (b) U/I and (c) V/I as a 
function of the length parameter / for varying direction of background mag- 
netic field; the magnitude \B\ and uo p are constant. Curves generated by di- 
rect numerical integration (solid) and adiabatic analytic calculation(dashed). 
Parameters gB = 1.0, L = 100, m^/ujp = 0.1, angles £(0) = 7r/2 and 
f (L) = vr/2 - 0.3tt; initial polarization (Q/I = 0, U/I = 1.0, V/I = 0.0). 
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2.1 Transition Probabilities 



Analytic calculations in the adiabatic limit fail for small frequencies, 
since in this case the transition probabilities between instantaneous 
eigenstates are large. Even in the large frequency regime the adia- 
batic limit fails unless the product gBL >> 1. This can be verified 
explicitly by computing the transition probabilities using the proce- 
dure discussed in Ref . JS] • The general solution to the wave equation 
can be written as 

\^ >= ^a n (z)e^o dz '^\ n > (45) 

n 

where |n > and to n are the instantaneous eigenmodes and eigenfre- 
quencies respectively. The evolution of the coefficients a n (z) with z 
gives an estimate of the transition among different eigenmodes. These 
coefficients are obtained by solving the equation 

d z b m ~ £ ^5*^^;^'^), (46) 

where we have approximated is defined by the equation 

a m {z) = e-y° dz ' (9 ^^b m (z) « b m (z) (47) 

In the small frequency regime the large transition probabilities are 
easily understandable: the mass matrix in Eq. I^has two eigenvalues 
very close to one another. In the opposite limit of large frequencies we 
find that the exponent in Eq. I4fil is small. This is because the exponent 
is inversely proportional tow, as long as gBL is small. The coefficient 
in Eq. ^2 for m = 1 and n = 2, for example, in this case is found 
to be proportional to Integrating this equation then gives a non- 
negligible contribution to the transition probability between different 
eigenstates. In the limit of large gBL >> 1 we find that the exponent 
is again large and suppresses the transition probabilities. 

Results: In Fig. Owe show a sample of results obtained in the case 
of varying direction of background magnetic field from analytic cal- 
culation in the adiabatic limit as well as direct numerical integration. 
Here angle £(0) = tt/2 and £(L) = 7r/2 — 0.3-7T, i.e. the transverse com- 
ponent of the background magnetic field is aligned along the y-axis 
initially and evolves to angle 0.3-7T after a distance L. The parameters 
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used in this figure are gB = 1.0, L = 100, m^/ujp = 0.1. The initial 
state of polarization has been chosen such that Q/I = 0, U/I = 1, and 
V/I = 0. The analytic results in this case are in good agreement with 
the numerical results, except in the limit of small frequencies. In the 
large frequency limit the exponent in Eq. \^is approximately equal 
to igBL/2. For the parameters chosen this phase factor is large and 
hence suppresses the transition probability between different eigen- 
states. 

In Fig. |3] we show a sample of results obtained in the case of 
varying direction of background magnetic field for a smaller value of 
the product gBL. Here we choose gB = 0.1, L = 100, and m^/u^ = 
0.1. In this case we use direct numerical integration since the analytic 
results are not reliable. The orientation of the background magnetic 
field is chosen to be same as in Fig. [21 i.e. £(0) = ir/2 and £(L) = 
7r/2 — 0.37T. The initial state of polarization has been chosen such that 
Q/I = 0, U/I = 1, and V/I = 0. The results obtained using this 
parameter choice and with uniform magnetic field direction are also 
shown for comparison. We find that the results obtained for the case of 
varying background magnetic field direction are considerably different 
in comparison to what is obtained in the case uniform direction. As 
expected the results agree in the limit of small uj. In Fig. |l]we show 
the results for the same parameter choice used in Fig. [21 but with the 
wave assumed to be unpolarized initially. 

The degree of polarization and the normalized Stokes parameters 
as a function of distance are shown in Fig. [3 Here the parameters 
are taken to be same as for the Fig. |3] with the length parameter 
I = 2u/(u)l - ml) = 10 and the wave is assumed to be unpolarized at 
source. We see that all the parameters p, Q/I, U/I, V/I oscillate with 
propagation distance. 

In Fig. El we show the relationship between Q/I and V/I for sev- 
eral different choice of parameters for the case of varying background 
magnetic field. The dependence of Q/I and V/I follows approximately 
an elliptical behaviour. This is in contrast to the the case of uniform 
magnetic field direction, which shows such a relationship between U/I 
and V/I |15| . As in the case of uniform background, a simple correla- 
tion is seen only for frequencies larger than a minimum frequency. At 
low frequencies the relationship becomes very complicated. 

One may be able to test relationships among different Stokes pa- 
rameters in future observations. To rule out other possible mecha- 
nisms affecting data, certain tests require observations over a suffi- 
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1000 




1000 



Figure 3: (a) Normalized Stokes parameters Q/I, U/I and V/I as a function 
of the length parameter I for varying background magnetic field direction; 
the magnitude \B\ and u p are constant. Parameters gB = 0.1, L = 100, 
rri^/Up = 0.1; angles £(0) = ir/2, = ix/2 - 0.37r; initial state of the 

polarization (Q = 0, U = 1.0, V = 0.0). Results for uniform background 
magnetic field (b) are shown for comparison. 




Figure 4: The degree of polarization p and the normalized Stokes parameters 
Q/I,U/I and V/I as a function of the length parameter I for varying direction 
of background magnetic field; the magnitude \B\ and u p are constant. Param- 
eters gB = 0.1, L = 100, m}/wj = 0.1; angles f (0) = tt/2, f (L) = 7r/2-0.37r. 
The wave is assumed to be unpolarized (Q — 0, /7 = V = 0) at source. 
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DISTANCE 



Figure 5: The degree of polarization p and the normalized Stokes parameters 
Q/I,U/I and V/I as a function of the distance of propagation for varying 
direction of background magnetic field; the magnitude \B\ and u p are con- 
stant. The parameters gB = 0.1, m^/u^ = 0.1, I = lujjiuj^ — mfy = 10; 
angles £(0) = 7r/2, £(L) = 7r/2 — 0.37T. The wave is assumed to be unpolar- 
ized (Q — 0, C/ = V = 0) at source. 
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Figure 6: A sample of results showing the correlation between the normalized 
Stokes parameters Q/I and V/I for some randomly chosen parameters and 
initial state of polarization. The results are shown for varying background 
magnetic field direction with the plasma frequency and the magnitude of 
the magnetic field uniform. Parameters (in arbitrary units) are (a) gB = 
2, L = 10, 0.04 < / < 20, (b) gB = 10, L — 10, 0.4 < I < 800, (c) 
gB = 1, L — 50, 0.2 <l < 800 and (d) gB = 10, L — 10, 0.04 < I < 100. 
The ratio mj/wj = 0.1; angles f (0) = vr/2 and £(L) = tt/2 - 0.3tt for all the 
plots. 
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ciently large frequency interval. 



3 Summary and Conclusion 

The general treatment of mixing of electromagnetic waves with pseu- 
doscalars in the presence of background magnetic field is a surprisingly 
intricate topic. The pseudoscalar mixes with (and indeed becomes) the 
longitudinal mode of light, a situation potentially generating cumula- 
tive deviation compared to treatments assuming the fields stay trans- 
verse. Cumulative errors do occur in principle, but for parameters of 
current interest they are fortunately controlled. The contribution due 
to the longitudinal component can be accommodated by redefining the 
pseudoscalar mass parameter — > m^ + ^Hi/e. This simplification 
led to exploring the problem of propagation in a magnetic field whose 
direction may vary along the path. The condition of adiabaticity is 
found to be rather stringent: For a wide range of parameter space the 
evolution cannot be assumed to be adiabatic. 

Thus the general problem of mixing of light with pseudoscalars has 
more twists and turns than could have been anticipated early. Stokes 
parameters show interesting correlations with one another which are 
distinctively different from those observed for fixed background field 
direction |15j. Such polarization effects may be observable with current 
technology, and may eventually serve either to identify new physics, 
or to put new limits on the pseudoscalar-photon coupling parameters. 

Acknowledgments: Work supported in part under Department of 
Energy grant number DE-FG02-04ER41308. 
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